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ABSTRACT 

We use a numerical simulation of cosmological reionization to estimate the likelihood of de- 
tecting Lyman-a emitting galaxies during the pre-reionization era. We show that it is possible to 
find galaxies even at z ~ 9 that are barely affected by the dumping wing of the Gunn-Peterson 
absorption from the neutral IGM outside of their H II regions. The damping wing becomes 
rapidly more significant at z > 9, but even at z > 10 is it not inconceivable (although quite hard) 
to see a Lyman-a emission line from a star-forming galaxy. 



Subject headings: cosmology: theory 
- galaxies: intergalactic medium 

Introduction 



cosmology: large-scale structure of universe - galaxies: formation 



A recent landslide in the number of detections 
of high redshift Lyman-a emitters (Dey et al. 1998; 
Hu, Cowie, & McMahon 1998; Weymann et al. 
1998; Hu, McMahon, & Cowie 1999; Stern & Spin- 
rad, 1999; Ellis et al. 2001; Ajiki et al. 2002; Hu 
et al. 2002; Kodaha et al. 2003; Taniguchi et al. 
2003; Rhoads et al. 2003; Maier et al. 2003; Cuby 
et al. 2003; Bunker et al. 2003; Hu et al. 2004; 
Pello et al. 2004; Kneib et al. 2004) offers a unique 
probe of the ionization history of the universe dur- 
ing cosmological reionization. As observations of 
Lyman-a absorption in the spectra of high redshift 
quasars (Djorgovski et al. 2001; Becker et al. 2001; 
Fan et a. 2002, 2003; White et al. 2003; Songaila 
2004) unambiguously indicate that the universe 
was reionized shortly before z — 6 (Gnedin 2002), 
the observations of Lyman-a at z > 6 directly 
probe the pre-reionization era. 

Specifically, since the damping wing of the 
Gunn-Peterson absorption (Gunn & Peterson 
1965) of the neutral IGM will affect the Lyman-a 



^CASA, University of Colorado, Boulder, CO 80309, 
USA; gnedin@casa.colorado.edu 

Ramon y Cajal Fellow, Institute de Astrofisica de An- 
dalucia (CSIC), E-18008 Granada, Spain; fprada@iaa.es 

^Visitor at the Kavli Institute for Cosmological Physics, 
University of Chicago 



emission line if the H II region around the emitting 
galaxy is not large enough (Miralda-Escude 1998; 
Haiman 2002), the observations of high redshift 
Lyman-a emitters tell us something about the dis- 
tribution and sizes of cosmological H II regions - 
although, to understand what specific constraints 
the observations place is a horrendously difficult 
task^ . 

However, the very fact of discovering a high- 
redshift Lyman-a emitter may potentially be used 
to constraint the ionization history of the universe. 
For example, Pello et al. (2004) report a tenta- 
tive discovery of a z = 10 galaxy. Is it possi- 
ble to see a Lyman-a emitting galaxy at such a 
high redshift? Loeb et al. (2004) used this fact 
to put some preliminary constraints on the ioniza- 
tion state of the IGM around that galaxy. While 
a real situation is likely to be more complicated 
because of the variety of reasons (H II regions are 
typically highly asymmetric because I-fronts prop- 
agate much faster across underdense voids then 
across dense filaments; bright galaxies are biased, 
so it is likely that more than one galaxy is lo- 
cated inside a single H II region; Lyman-a emit- 
ters can also be located inside H II regions of lumi- 



^ Since most Lyman-a emitters are discovered because they 
are highly lensed by foreground clusters of galaxies, de- 
termining the observational selection function is extremely 
difficult. 
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nous quasars, which arc often many times larger 
than H II regions of galaxies, etc) , it is interest- 
ing to investigate how likely it is to find a high 
rcdshift galaxy unaffected by the damping wing of 
the Gunn-Peterson absorption from the surround- 
ing IGM. 

2. Star Formation Modeling 

In order to construct a model for the Gunn- 
Peterson absorption in the pre-reionization era, we 
use a numerical simulation of cosmological reion- 
ization similar to the one described in Gnedin 
(2000). The new simulation, however, incor- 
porates a significantly more accurate method 
for following the time-dependent and spatially- 
inhomogeneous radiative transfer using the Opti- 
cally Thin Variable Eddington Tensor (OTVET) 
approximation of Gnedin & Abel (2001). 

The simulation assumes a flat cosmology with 
the values of cosmological parameters as measured 
by the WMAP satellite (Spcrgcl ct al. 2003)^. The 
size of the simulation volume is 8h~^ Mpc in co- 
moving units, and the comoving spatial resolution 
of the simulation is 2h~^]ipc. The simulation in- 
cluded 128^ dark matter particles (with mass of 
2.6 X IO^Mq), an equal number of baryonic cells 
on a quasi-Lagrangian moving mesh, and about 3 
million stellar particles that formed continuously 
during the simulation. Only PopII stars are in- 
cluded in the simulation as sources of ionizing ra- 
diation. 

The simulation was adjusted to reproduce the 

observed evolution of the mean transmitted flux 
in the Lyman-a line between z = 5 and z = 6.4 
(White et al. 2003; Songaila 2004). We, therefore, 
can be confident that the evolution of the IGM 
is treated correctly in the simulation at least at 
~ 6. 

Given a simulation, we identify all star-forming 
galaxies using the DENMAX halo finding algo- 
rithm of Bertschinger & Gelb (1991). For each 
galaxy, we generate 100 lines of sight that origi- 
nate at the galaxy location and go along randomly 
chosen directions. We then generate an absorp- 
tion spectrum along each line of sight following 
the standard procedure. For the computational 

2 Specifically, we assume Qm = 0.27, f2j, = 0.04, h = 0.71, 
n = 1, and as = 0.85. 
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Fig. 1. — Baryon number density (black lines) 
and neutral hydrogen number density (gray lines) 
along 3 random lines of sight originating at the 
location of the most luminous galaxy at z = 9 in 
our simulation. 

efficiency, we use Liu et al. (2001) approximation 
for the Voigt profile. 

Figure 1 illustrates the distribution of the 
baryon number density and neutral hydrogen 
number density along 3 random lines of sight origi- 
nating at the location of the most luminous galaxy 
at z — 9. The H II region around the central galax- 
ies is easily visible. Because the brightest galaxies 
are highly clustered, smaller galaxies and other 
dense structures are often located within the H II 
regions of brighter ones. For example, in the top 
panel, an overdensity of the order of 10 is located 
within the H II region of the central galaxy. This 
overdensity is the outer part of a smaller galaxy, 
and is, therefore, significantly more neutral than 
the rest of the H II region (in fact, more neu- 
tral than simple ionization equilibrium would pre- 
dict, due to shadowing and time-dependent ion- 
ization). This neutral spike is sufficient to produce 
large damping wing absorption at the location of 
the galaxy and make the Lyman-a of the central 
galaxy unobservable. Two bottom panels, how- 
ever, show different, "clean" lines of sight that are 
not obscured by neighbors and produce lower ab- 
sorption by the damping wing of the neutral IGM 
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outside the H II region. This iUustratcs that the 
effect of the damping wing of the Gunn-Peterson 
absorption strongly depends on the direction of 
the hne of sight: the same Lyman-a emitting 
galaxy may be observable or non-observable de- 
pending on the direction of viewing. 

A significant complication, however, for the 
evaluation of the strength of the Lyman-a emis- 
sion line from the central galaxy is the galaxy it- 
self. Synthetic Lyman-a spectra from simulation 
are often treated as absorption spectra, when, in 
reality, they are scattering spectra. For gas along a 
line of sight well outside the source of emission this 
is an appropriate approximation, since scattering 
off the direction of viewing removes the photons 
that could be observed, and, thus, appears as effec- 
tive absorption. This is not true for the source of 
emission, since photons that were originally emit- 
ted in directions different from the line of sight 
are scattered into the direction of viewing. Thus, 
without the proper line radiative transfer treat- 
ment (which is not included in our simulation), 
it is not possible to accurately predict shapes of 
the Lyman-a emission lines from simulated galax- 
ies. Leaving such treatment for future work, in 
this paper we are concerned only with an-order- 
of- magnitude estimate of the effect that damping 
wings of the neutral gas outside H II regions have 
on Lyman-a emission lines. 

In order to obtain such an estimate, we adopt 
a following simple procedure. We assume that the 
Lyman-a emission line from a typical galaxy has 
a width of the order of 150 — 200km/s. Thus, we 
exclude all Lyman-a absorption within a specified 
velocity distance Woff from the systemic velocity of 
a simulated galaxy (because within this distance 
we should treat the Lyman-a as scattering and not 
as absorption). Also, since the significant suppres- 
sion of Lyman-a emission at the systemic velocity 
may still make the red wing of the emission line 
visible, we choose the Lyman-a optical depth at 
— 150km/s as a quantity that roughly character- 
izes the damping wing absorption of the Lyman- 
a emission line. Typical Lyman-a emission lines 
have gaussian widths of the order of 60 — 80 km/s, 
so — 150km/s roughly corresponds to 2a to the 
red of the line systemic velocity. In this paper we 
use a range of values for VoS to estimate the sen- 
sitivity of our results to that somewhat arbitrary 
parameter. 
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Fig. 2. — The limiting Lyman-a fiux decrement as 
a function of star formation rate of a central galaxy 
for top 3% (solid lines) and top 10% (long-dashed 
lines) of all galaxies. Black lines show the case 
of VoS = 150 km/s, and gray lines show the case 
of Uofi' = 200 km/s. The short-dashed and dotted 
lines show the the limiting Lyman-a flux decre- 
ment as a function of star formation rate for iso- 
lated galaxies with 100 kpc and 300 kpc (in phys- 
ical units) exclusion radius respectively (see text 
for more details). 

With this approach, we can calculate the 
Lyman-a optical depth at —150 km/s for each of 
100 lines of sight for each of the simulated galaxy. 
Thus, there are 100 values of the damping wing 
absorption optical depth TLya, or, equivalently, 
the flux decrement Fi,ya = exp(— TLya), for each 
galaxy. These 100 values are not all the same, 
of course, so we can use them to approximately 
characterize the probability that a given galax;y 
(labeled by its star formation rate) has a given 
value of the Lyman-a flux decrement. 

3. Results 

As a test of our approach, we show in Figure 2 
the values for the Lyman-a flux decrement above 
which 3% and 10% of all galaxies lie - for example, 
solid lines indicate that for 3% of all galaxies the 
decrease in flux at —150 km/s is not more than 



DRAFT: February 2, 2008 



4 



a factor of 10. The two values for the parameter 
VoS are used - the difference between black and 
gray lines roughly indicates the uncertainty of our 
estimate. 

An interesting feature of Fig. 2 is the lack of de- 
pendence of the flux decrement on the star forma- 
tion rate. Naively, one would expect that galax- 
ies with lower star formation rates create smaller 
H II regions around them, and, therefore, should 
exhibit higher suppression of their Lyman-a emis- 
sion lines - clearly at odds with our result. How- 
ever, this naive argument misses an important fact 
that high rcdshift galaxies arc highly clustered, so 
it is likely for a faint galaxy to be located close 
to a bright galaxy, within the bright galaxy H II 
region. In order to test that this effect is respon- 
sible for the lack of dependence of the flux decre- 
ment on the star formation rate, we compute the 
flux decrement versus the star formation rate for 
the sub-sample of isolated galaxies. We define a 
galaxy to be isolated if it does not contain a galaxy 
with a higher star formation rate within a sphere 
of predefined radius risoi- Short-dashed and dot- 
ted lines in in Fig. 2 show the flux decrement for 
the top 10% of all galaxies (an analog of the long- 
dashed line) for sub-samples of isolated galaxies 
with risoi = 100 kpc and 300 kpc (in physical units) 
respectively. As one can see, isolated galaxies ex- 
hibit a strong dependence on the star formation 
rate for rates above about O.lM0/yr. At even 
lower star formation rates a different effect takes 
over: since the star formation in dwarf galaxies 
is highly episodic (c.f. Ricotti, Gnedin, & ShuU 
2003), the instantaneous value of the star forma- 
tion rate and the size of an H II region cease to 
correlate. 

The redshift dependence of the Lyman-a flux 
decrement is shown in Figure 3. As once can see, 
our estimates become highly uncertain at z ^ 10 - 
at these redshifts full radiative transfer in Lyman- 
a line is required to make an accurate prediction 
for the effect of the damping wing on the profile 
of the Lyman-a emission line. At lower redshifts, 
however, our estimates are more precise. 

4. Discussion 

Somewhat unexpectedly, even at as high a red- 
shift as z ^ 9, one percent of all Lyman-a emitters 
should have their emission lines barely affected by 




Fig. 3. — The limiting Lyman-a flux decrement 
as a function of star formation rate of a cen- 
tral galaxy for top 1% (dotted lines), 3% (dashed 
lines), and top 10% (solid lines) of all galaxies as 
a function of redshift. Gray lines show the case of 
t^off = 200km/s, and black lines show the case of 
= 150km/s - the difference between the cor- 
responding lines indicates the uncertainty of our 
calculations. 

the damping wing. At higher redshifts the damp- 
ing wing quickly becomes more significant, but 
finding Lyman-a emitters at z > 10 should still 
be possible, although they must be quite rare. 

The reason for this rapid transition is quite sim- 
ple: as the universe expands, H II regions increase 
in size, so at each moment there exists a character- 
istic scale for the size distribution of H II regions 
(say, a typical size of an H II region around an 
galaxy). At the same time, the damping wing of 
the Gunn-Peterson absorption introduces another 
spatial scale (about Az = 0.01 or proper distance 
of 200 kpc at z ~ 10; Miralda-Escide 1998). Thus, 
when typical H II regions around bright galax- 
ies exceed this scale, a substantial fraction of all 
Lyman-a emitters becomes observable. 

It is important to emphasize here that our es- 
timate is likely only a lower limit: because of the 
finite size of the computational box, we miss the 
brightest galaxies which create the largest H II re- 
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gions (although the size of the simulation box was 
selected so as to make this effect not significant 
at ^; > 6). In addition, our simulation does not 
include quasars, which can create H II regions sev- 
eral times that of galaxies. High redshift quasars 
are rare (there should be only about 2 SLOAN 
quasars in the whole observable universe between 
redshifts 9 and 10, Fan et al. 2001), so that effect 
is most likely not large too, but without a detailed 
simulation we cannot estimate it, so we treat our 
results as a lower limit. 

We thank A. Kravtsov for valuable comments 
and the Kavli Institute for Cosmological Physics 
at the University of Chicago for hospitality dur- 
ing the time this paper was initiated. This work 
was supported in part by NSF grant AST-0134373 
and by National Computational Science Alliance 
under grant AST-020018N and utilized IBM P690 
array at the National Center for Supercomputing 
Applications. 
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